Formation of cosmological mass condensation within a FRW 
universe: exact general relativistic solutions 
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Abstract 

Within the framework of an exact general relativistic formulation of gluing mani- 
folds, we consider the problem of matching an inhomogeneous overdense region to a 
Friedmann-Robertson- Walker background universe in the general spherical symmet- 
ric case of pressure- free models. It is shown that, in general, the matching is only 
possible through a thin shell, a fact ignored in the literature. In addition to this, in 
subhorizon cases where the matching is possible, an intermediate underdense region 
will necessarily arise. 
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Large scale structure formation is a challenging field from both theoretical and 
observational points of view. Motivated by the observational data on the existence of 
large scale voids in the universe, Sato and co-workers considered in a series of papers 
an underdense spherical region immersed in a Friedmann-Robertson- Walker (FRW) 
background universe applying the thin shell formalism of general relativity [1] . On 
the other hand, the appearance of holes with densities less than the average cos- 
mological density around rich spherical clusters of galaxies has been seen by some 
authors [2-4] . In this paper, we model such inhomogeneous large scale structures as 
a Lemaitre-Tolman-Bondi (LTB) manifold with different density profiles glued to a 
homogeneous pressure-free FRW background from which a sphere of dust matter is 
removed. Our calculation is based on an exact general relativistic formulation of 
gluing manifolds. We may therefore consider our contribution as a generalization 
of the work done by Olson and Silk[3]. Their calculation was made basically within 
Newtonian dynamics without being cautious about the matching conditions. 
Consider first a spherical inhomogeneous region containing dust matter, represented 
by a LTB cosmological model embedded in a pressure-free FRW background uni- 
verse with the uniform density pb- We choose the LTB metric to be written in the 
synchronized comoving coordinates in the form [5] 

ds 2 = -dt 2 + — t dr 2 + R 2 (r,t)(d6 2 + sin 2 9d<J) 2 ). (1) 
1 + E(r) 

The overdot and prime denote partial differentiation with respect to t and r, respect- 
ively, and E(r) is an arbitrary real function such that E(r) > —1. Then the corre- 
sponding Einstein field equations turn out to be 

R\r,t) = E(r) + 2 -^f±, (2) 
M'(r) 

To be as general as possible, we let the density p(r, t) to be a general function of r, 
i.e., it is not necessarily a monotonic function of radial distance from the center of 
the region with mass condensation. M(r) is defined by 

Mir) = 4tt p(r, t)R 2 dR- (4) 

Jo 

Furthermore, in order to avoid shell crossing of dust matter during their radial 
motion, we must have R'(r, t) > 0. The solution to the above equations shows that 
an overdense spherical inhomogeneity with E(r) < within R evolves just like a 
closed universe, namely it reaches to a maximum radius at a certain time, then 
the expansion ceases and undergoes a gravitational collapse so that a bound object 
forms in such a way. 

Now, in the literature on large scale structures of the universe it is always tacitly 
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assumed that the junction of an overdense region to a background FRW universe is 
a continuous one and there is no need of a singular hypersurface along which the 
gluing is made[6]. Let us denote by E the (2+l)-dimensional timelike boundary of 
the two distinct spherically symmetric regions glued together. We will show that 
this is in general not true and in the most cases of interest E cannot be just a 
boundary surface but is a singular hypersurface carrying energy momentum. 
We now write down the appropriate junction equation on E expressing the jump of 
the angular component of extrinsic curvature tensor K$ across E as [7,8] 





^iP I AnaR , (5) 



where = means that both sides of the equality are evaluated on E, r is the proper 
time of the comoving observer on E, and a is the surface energy density of the shell 
E. The sign functions are fixed according to the convention ei n (e out ) = +1 for R 
increasing in the outward normal direction to E, while e in {e out ) = — 1 for decreasing 
R. An average density for the inhomogeneous region is also defined as 



P = 



M(r) 



f R 3 



(6) 



Note that in general the matching is only possible if a thin layer is formed on the 
boundary of the two manifolds, i.e. a ^ 0, except for the case where p = pb, and 

(in ^out- 

Now, similar to the approach used by Berezin et al [8], we solve Eq. (5) for 
to obtain 



(7) 



where £ is defined by 



^ < 8 > 



Noting that the surface energy density a on E is positive, we substitute Eq. (7) 
back into Eq. (5) to get 

Q n |C + l| -e out |£-l| = 2- (9) 

It is easily seen that the sign functions e in / out are determined by the value of £: 

:+l,+l) for £>1, 
(e in ,e out ) = { (+1,-1) for |f| < 1, (10) 
(-1,-1) for £<-!• 
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Now, e out may be related to different parameters of FRW universe. It is known 
that [9] 



tout = sgn h v b H b R 




(11) 



where H b is the Hubble parameter of the background, v b is the peculiar velocity of 
E relative to the background. The two coordinates r b and \b are related to each 
other as 



Depending on the density distribution within the shell, the following cases may now 
be distinguished: 

(I) Junction of an overdense region to an background FRW universe, i.e., p(r, r) > p b . 
Obviously £ < 0. Therefore, according to Eq. (10), we must have e out = — 1. This 
leads to the following cases: 

(1) k = 0, —1; R < H^ 1 . In both these cases it is easily seen from Eqs. (11), (12) 
that e out — +1. Therefore, there is no matching of an overdense region smaller than 
the Hubble radius to a background FRW universe with k — 0, — 1. 

(2) k — 0, — 1; R > H^ 1 , and the peculiar velocity v b < 0. It is seen from Eq. (11) 
that e out = — 1[9]. Therefore, overdense structures with a radius greater than the 
Hubble radius and v b < may be glued to a background FRW with k — 0, — 1. It is 
trivially seen that there is no matching for the case of v b > 0. 

(3) k = +1. Taking the definition Eq. (12), we see from (11) that it is possible to 
have e out — — 1. Therefore the matching is possible without any restriction. 
Except from this last case of a closed FRW background universe and the case of a 
superhorizon overdense region discussed in (2), we have seen that it is not possible to 
match an overdense spherical region to a background FRW universe. Now, to look 
for more possible matching let us assume the overdense structure to be surrounded 
by an underdense region. Three other cases may be distinguished now. 

(II) p < p b ; k — 0, — 1. We have therefore an overdense inhomogeneous region sur- 
rounded continuously by an underdense region p(r,t) < p b , i.e., a void described by 
the same LTB metric. Now, note that the junction of the LTB region at the under- 
dense void part of it to the background FRW is also described by Eq. (5), while p is 
the overall average density related to the total mass enclosed within the boundary 
separating the hole region of physical radius R from the uniform background. 

(1) R < H b ~ l . It is seen from Eq. (11) that e out = +1, irrespective of the sign of the 
peculiar velocity. Therefore, according to Eq. (10), the matching is possible with 
£ > 1 via a thin shell having a ^ 0. 

(2) R > H b _1 . In this case for the negative peculiar velocity we have e out = — 1. 
Therefore the matching is possible with < £ < 1, i.e., a ^ 0. If the peculiar 
velocity is positive we are faced with a case exactly the same as (1). 




sinxb (k = +1, closed universe) 
Xb (k = 0, flat universe), 
sinh^b (k — — 1, open universe)- 



(12) 
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(III) p = pb, k = 0, — 1. Now, the void completely compensates the overdense region, 
so that the overall mean density is equal to the background density. 

(1) R < H b ~ x . It turns out that a = 0, i.e., the matching is possible without forma- 
tion of a thin shell. This situation reminds us of the Einstein-Strauss type model 
[10] where the overdense region is surrounded by the vacuum shell to compensate 
the mass excess. 

(2) R > H b l . For the negative peculiar velocity, we have again e out = —1. Both 
junctions with a = and ^ are possible, as it is easily seen from Eq. (5). The 
case €i n = — 1 corresponds to a = meaning no thin shell, and the case €i n = +1 is 
only possible through a thin shell. The case of positive peculiar velocity is similar 
to the (1) case. 

(IV) p > pi,; k — 0, — 1. This is similar to the matching where no void is in between, 
as discussed in the 

cases considered in (I). 

Table I summarizes the above results. 

Table I. Classification of junctions of LTB spherical inhomogeneity to the fiat and open 
FRW background. An extensive version of the paper with the detailed calculation and 
critical review of the relevant literature is going to be published in [11]. 





R > H b ~\ v p <0 

{tout — — 1) 


R < H b ~ L or 
R > H b -\ v p >0 

(tout = +1) 


p> Pb 


Thin shell 


Junction impossible 


p = pb 


No thin shell 
(l) e --1 


No thin shell 

tin = +1 


/..\ Thin shell 

<"> e = o 


p< pb 


Thin shell 

o<e<i 


Thin shell 
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